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annexin Al; LGALS1, beta-galactoside-binding lectin; PRDX, peroxiredoxin; NUMA, 
nuclear mitotic apparatus protein; LMNA, lamin-A; EIF5A, eukaryotic translation initiation 
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PSMB2, proteasomel subunit beta type-2; ANXA2, annexin A2; HSPB1, Heat shock protein 
27; GST, glutathione S-transferase; SOD1, Cu/Zn superoxide dismutase. 
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Abstract 

Infectious bronchitis coronavirus (IBV), Newcastle disease virus (NDV), and avian 
influenza virus (AIV) H9 subtype are major pathogens of chickens causing serious 
respiratory tract disease and heavy economic losses. To better understand the replication 
features of these viruses in their target organs and molecular pathogenesis of these different 
viruses, comparative proteomic analysis was performed to investigate the proteome changes 
of primary target organ during IBV, NDV and AIV H9 infections, using 2D-DIGE followed 
MALDI-TOF/TOF-MS. In total, 44, 39,41, 48 and 38 proteins were identified in the tracheal 
tissues of the chickens inoculated with IBV (ck/CH/FDF/97I, H120), NDV (Fa Sota), and 
AIV H9, and between ck/CH/FDF/97I and HI20, respectively. Bioinformatics analysis 
showed that IBV, NDV and AIV H9 induced similar core host responses involved in 
biosynthetic, catabolic, metabolic, signal transduction, transport, cytoskeleton organization, 
macromolecular complex assembly, cell death, response to stress and immune system process. 
Comparative analysis of host response induced by different viruses indicated differences in 
protein expression changes induced by IBV, NDV, and AIV H9 may be responsible for the 
specific pathogenesis of these different viruses. Our result reveals specific host response to 
IBV, NDV, and AIVH9 infections and provides insights into the distinct pathogenic 
mechanisms of these avian respiratory viruses. 
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1 Introduction 

Avian respiratory viruses, including infectious bronchitis coronavirus (IBV), Newcastle 
disease virus (NDV), and avian influenza virus (AIV), are the primary cause of morbidity and 
mortality worldwide in the poultry industry and have been identified as the most 
economically important etiological agents of acute, highly contagious poultry diseases. 
Among these viruses, NDV is a member of the Avulavirus genus within the Paramyxoviridae 
family and is prevalent worldwide and infects almost all avian species. In fact, highly 
contagious respiratory, enteric, and neurological diseases in chickens are often caused by a 
variety of infectious NDV strains [1,2]. AIV is an enveloped virus with a segmented, single- 
stranded, negative-sense RNA genome belonging to the Orthomyxoviridae family and can 
cause clinical diseases in a wide range of mammalian and avian species. Low pathogenic 
avian influenza virus subtype H9 circulates worldwide in multiple avian species thereby 
posing a potential threat to not only the poultry industry, but also to public health [3,4]. IBV 
is a gammacoronavirus of the subfamily Coronavirinae, family Coronaviridae, and order 
Nidovirales [5], and primarily replicates in epithelial cells of the respiratory tract and causes 
lower-respiratory tract infections in chickens ( Gallus domesticus). In addition, emerging IBV 
isolates that vary in their replicative abilities in non-respiratory tissues have been associated 
with proventriculus and kidney lesions [6, 7]. 

Proteomic analyses of pathogen-infected organisms can offer significant insight to the 
“battle” between pathogens and host during infection, hence proteomic analyses of host 
cellular responses against viral infections could aid in identification of potential host factors 
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involved in the viral life cycle and host defense mechanisms. To date, various proteomic 
strategies have been applied to elucidate protein dynamics during influenza infections, such 
as differential expression of brain proteins in response to the neurovirulent AIV H5N1 
infection [8], host responses of murine lung tissues infected with AIV H5N1 [9], and cellular 
protein changes in H1N1 -infected human lung A549 cells [10]; while cellular responses to 
AIV H9N2 subtype infections has only been reported by Liu et al. in human cells [11]. For 
IBV, proteomic changes have been reported in the cytoplasm, nuclei, and nucleoli in 
mammalian Vero cells and nucleolar proteomes in avian cells (DF-1) infected with IBV 
(Beaudette strain) [12, 13]. We previously reported alterations in the proteome of tracheal and 
kidney cells during IBV strain HI 20 i nf ection in vitro, and IBV strain ck/CH/LDL/97I and its 
attenuated strain ck/CH/LDL/97I Pus infection in vivo [14, 15]. However, alterations of 
cellular protein expression profiles, whether in NDV-infected cell lines or tissues, and the 
differences in induced protein expression changes between different avian respiratory viruses 
(IBV, NDV and AIV H9) have not yet been reported. 

The tracheal epithelium of chicken is not only the primary target for replication and 
infection of avian respiratory viruses, but also the first line of defense in the innate immune 
response. The avian respiratory viruses, including IBV, NDV and AIV H9, can infect tracheal 
epithelial cells and induce damage to tracheal epithelium resulting in respiratory disease. 
However, the severity and outcome of the clinical illness caused by these different viruses are 
different [16-18]. To better understand the molecular basis of pathogenesis and specific virus- 
host interactions, herein, we described the alterations of protein expression profiles in chicken 
tracheal tissues caused by IBV ck/CH/LDL/97I and HI 20, NDV La Sota, and AIV H9 
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infection by using 2D-DIGE analyses. A series of differentially expressed proteins associated 
with viral life cycle and host response to IBV, NDV and AIV H9 infections were identified. 
The present proteome results using the same virus of different strains (IBV virulent 
ck/CH/LDL/97I and HI20 vaccine) and different viruses (IBV, NDV and AIV H9) might 
provide more significant information on specific aspects of molecular pathogenesis and virus- 
host interactions of these different viruses. Furthermore, these findings might also provide 
useful clues for the development of novel prevention or therapeutics strategies against these 
viruses. 
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2 Materials and methods 

2.1 Viral strains 

A commercial vaccine against IBV strain HI 20 was obtained in freeze-dried 1000-dose 
vials that each contained at least 10 3 egg infective dose 50% (EID 50 ) per dose (isbi® BIO 
H120; batch No. F48355). IBV strain ck/CH/LDL/97I [19], NDV strain La Sota vaccine, and 
AIV subtype H9 [20] were used in this study. The virus stock for this study was produced by 
inoculating the virus into embryonated SPF chicken eggs via the allantoic cavity and 
collecting the infectious allantoic fluid 72 h post-inoculation. The allantoic fluid was clarified 
by centrifugation at 3000 x g for 10 min and filtered with a Teflon membrane and the 
presence of IBV was confirmed by negative contrast electron microscopy (JEM-1200, EX; 
JEOL Ltd., Tokyo, Japan) and reverse transcription polymerase chain reaction (RT-PCR) as 
previously described [21]. NDV and subtype H9 AIV viral particles in the allantoic fluids of 
inoculated eggs were confirmed using hemagglutination inhibition (HI) analysis with specific 
antibodies, respectively [22], The titers of the viruses were determined by inoculation at 10- 
fold dilutions into groups of five 10-day-old embryonated chicken eggs. The EID 50 was 
calculated by the method of Reed and Muench [23]. 
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2.2 Experimental design 

One hundred and twenty 1-day-old SPF White Leghorn chicks were divided into five 
groups of 24 birds each and housed in different isolators with negative pressure. Chickens 
assigned to groups 1-4 were inoculated with IBV ck/CH/LDL/97I, H120, NDV La Sota, and 
AIV H9, respectively, by oculonasal application at 3 days of age with a dose of loglO 5 
EID 50 /IOO pLper chick. Birds in group 5 were mock-inoculated with sterile allantoic fluid as 
a negative control. This study was approved by the Animal Welfare Committee of 
Heilongjiang Province, China. All experiments were conducted at appropriate biosafety levels. 

Three birds from each group were killed humanely at 3, 7, 12, and 21 days post¬ 
inoculation (dpi) and fresh tissue of trachea samples were quickly harvested and washed with 
ice-cold Tris-sucrose buffer (10 mM Tris, 250 mM sucrose, pH 7.5), snap-frozen in liquid 
nitrogen, and stored at-80°C for subsequent 2D-DIGE analysis. Sera samples from 10 
chickens in each group were collected at 3, 7, 12, and 21 dpi for enzyme-linked 
immunosorbent assay (ELISA) and HI analyses. Sera samples from groups 1,2, and 5 were 
assayed using a commercial IBV antibody test kit (IDEXX Corporation, Westbrook, Maine, 
USA) according to the manufacturer’s instructions. Semm samples from groups 3, 4, and 5 
were detected for HI antibody titer (positive when HI >3) against NDV and AIV, respectively, 
using an HI assay as reported previously [24]. Viral load of IBV, NDV and AIV in tracheal 
tissues were analyzed using methods described previously. Briefly, Viral RNA was extracted 
from the tracheal tissues of sacrificial birds from each group using TRIzol Reagent 
(Invitrogen, Carlsbad, CA, USA). For detection of the viruses, specific primers and TaqMan 
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probes against IBV [25], NDV [26], and AIV subtype H9 [4] viral RNA were selected. Real¬ 
time RT-PCR was performed on a LightCycler® 480 real-time PCR system with One Step 
PrimeScript™ RT-PCR Kit (TaKaRa, Dalian, China) [27], All of the samples were tested in 
triplicate in each reaction. Standard templates with known concentration and no template 
negative control were contained for quantitative analysis. 

2.3 2D-DIGE analysis of tracheal proteins 

The tracheal tissues were ground into powder in liquid nitrogen. Fifty milligrams of 
powder were extracted with 500 pL of DIGE-compatible lysis buffer containing 7 M urea, 2 
M thiourea, 4% CHAPS, 1% nuclease mix, and 1% Protease Inhibitor Mix (GE Healthcare, 
Waukesha, WI, USA), vortexed for 3 h in an ice bath, and then centrifuged at 15 OOOxg for 1 
h at 4°C. The supernatant was collected and purified using the Ettan 2-D Clean-Up Kit (GE 
Healthcare). The pH of the protein sample was adjusted to 8.5 and the protein concentration 
was determined using the Ettan 2-D Quant Kit (GE Healthcare). The purified protein lysates 
(50 pg) were labeled with 400 pmol of Cy2, Cy3, and Cy5 dyes (GE Healthcare) following 
the manufacturer’s instructions. Briefly, the Cy2-labeled internal standard was created by 
pooling an aliquot from all involved biological samples, while Cy3 and Cy5 were used 
randomly to label samples from the IBV ck/CH/LDL/97I-, H120-, NDV La Sota-, AIV H9-, 
and mock-infected groups. The experimental design of the different fluorescent dye labeling 
protocols is listed in Supporting Information Table S1. The labeling reaction was performed 
as described in our previous study [15] and labeled samples were subsequently loaded onto 
24-cm, linear pH 3-10 IPG strips (GE Healthcare) for integrated isoelectric focusing (IEF) 
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which was performed using a Ettan™ IPGphor 3 fully integrated IEF system (GE Healthcare). 
After IEF, the IPG strips were treated in an equilibration buffer (6 M urea, 50 mM Tris-HCl, 
pH 8.8, 30% glycerol, 2% SDS, and 0.002% bromophenol blue) containing 1% dithiothreitol 
and 2.5% iodoacetamide. A second dimension sodium dodecyl sulfate polyacrylamide gel 
electrophoresis (SDS-PAGE) was performed using an Ettan DALT six system electrophoresis 
system (GE Healthcare). 

After SDS-PAGE, the CyDye-labeled gels were scanned using the Typhoon 9400 
scanner (GE Healthcare) at a resolution of 100 pm. The intensity was adjusted to ensure that 
the maximum volume of each image was within 60,000-80,000. DIGE analysis was 
performed using DeCyder software version V6.5 (GE Healthcare). Protein spots were 
pairwise compared among control and different infection groups and differentially expressed 
spots with statistically significant results via the Student’s /-test (p < 0.05) and average ratio 
>1.2 or <-1.2 were selected for MS identification. Spot picking was performed using 
preparative gels under the conditions as described above for 2D-DIGE except that the IPG 
strips were loaded with 1000 pg of protein and the gels were stained with Coomassie Blue R- 
350, then the protein spots were excised manually from the gels and subjected to MALDI- 
TOF/TOF MS identification. 
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2.4 MALDI-TOF/TOF MS and database search 

The gel spots were excised, trypsin-digested, and subjected to MALDI-TOF/TOF MS 
analysis as described in our previous study [15]. The peptide mass fingerprint (PMF) 
combined MS/MS data were submitted to MASCOT version 3.0 (Matrix Science, Inc., 
Boston, MA, USA) for identification according to the National Center for Biotechnology 
Information (NCBI) nr database (release 16/01/2010, 10,343,571 sequences, and 
3,528,215,794 residues). The search parameters were set as follows: Gallus, trypsin cleavage 
(one missed cleavage allowed), carbamidomethylation of cysteine as a fixed modification, 
oxidation of methionine as a variable modification, peptide mass tolerance set at 100 ppm, 
and fragment tolerance set at 0.8 Da. The criterion for successful identification of a protein 
was a protein score confidence interval (Cl %) >95%. 

2.5 Gene ontology (GO) annotation and classification 

To investigate the potential roles played by the differentially expressed proteins during 
infection, GI numbers of identified proteins were submitted to GORetriever 
( http://www.agbase.msstate.edu/ ) to obtain the Gene Ontology (GO) annotations. If no 
annotation was acquired from GORetriever, GOanna was used to retrieve GO annotations 
assigned on the basis of searches using the Basic Local Alignment Search Tool. Then, the 
resulting annotations were summarized using GOSlimViewer. 
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2.6 Analysis of mRNA abundance by real-time RT-PCR 

Ten target proteins which extensively changed during IBV, NDV and AIV H9 infections 
shown by 2D-DIGE and potentially involved in the viral infection were selected for 
transcriptional analysis. The expression of selected proteins were not all altered at four time 
points (Supporting Information Table S2), therefore the mRNA levels of selected proteins at 
their corresponding time point were analyzed, respectively. Specific primers for the target 
genes of MS/MS-identified proteins were designed according to the published gene 
sequences using Beacon Designer software 7.5 (Premier Biosoft International, Palo Alto, CA, 
USA) and information on the primers is listed in Supporting Information Table S3. Total 
tracheal tissue RNA was extracted from sacrificed birds of each group at each time point 
using TRIzol Reagent (Invitrogen) and then reverse-transcribed to cDNA. The real-time PCR 
was performed using a LightCycler® 480 real-time PCR system (Roche Diagnostics 
Deutschland GmbH, Mannheim, Germany). Melting curve and quantitative analyses of the 
data was determined using the LightCycler® 480 real-time PCR system with the relative 
quantification 2~ AACt model, Cp values of target gene in control group at each time point were 
used as the calibrator (relative expression=l), and 28S ribosomal RNA gene was used as an 
internal reference gene. 
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2.7 Western blot analysis 

To verify the 2D-DIGE results, several representative proteins were selected for Western 
blot analysis based on the commercial availability of the corresponding antibodies. Tracheal 
proteins from sacrificed birds of each group at all four time points were extracted and 
equivalent lysate amounts were subjected to 12% SDS-PAGE and then transferred to 
nitrocellulose membranes. After blocking overnight at 4°C, the membranes were incubated 
with goat polyclonal antibody to annexin A2 (ANXA2; sc-1924; Santa Cruz Biotechnology), 
rabbit polyclonal antibody to annexin A1 (ANXA1; sc-11387; Santa Cruz Biotechnology), 
goat polyclonal antibody to peroxiredoxin I (PRDX1; sc-23970; Santa Cruz Biotechnology), 
mouse monoclonal antibody to heat shock protein 27 (HSPBl; sc-51956, Santa Cruz 
Biotechnology), rabbit polyclonal antibody to phospho-HSPBl (pSer82) (SAB4504434, 
Sigma-Aldrich) and mouse monoclonal antibody to [S-tubulin (sc-55529; Santa Cruz 
Biotechnology) for 1 h at 37°C, respectively. Then, the membranes were incubated with 
IRDye800DX 'conjugated anti-mouse, -goat, and -rabbit secondary antibodies (610-130-121, 
605-430-002, and 611-730-127; Rockland Immunochemicals, Gilbertsville, PA, USA), 
respectively, for 1 h at 37°C. After washing, the membranes were scanned using a LI-COR 
infrared imaging system with Odyssey software (Li-Cor Biosciences, Lincoln, NE, USA). 
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3 Results 

3.1 Clinical observations, serology and virus detection in trachea 

All chickens infected with the IBV ck/CH/LDL/97I strain exhibited clinical respiratory 
signs at -3-12 dpi and one of the diseased birds died at 12 dpi. The clinical signs included 
tracheal rales, watery eyes, nasal mucus, and sneezing. Gross lesions within the chickens 
were confined mainly to the kidneys [15], whereas the chickens in the control group and 
those inoculated with the IBV HI 20 and NDV La Sota strains showed no apparent respiratory 
clinical signs and no gross lesions. AIV H9-infected chickens appeared lethargy, depression, 
coughing, sneezing and dyspnea during the experimental period. 

The serological responses induced by infections with IBV ck/CH/LDL/97I, HI20, NDV 
La Sota, and AIV H9 are summarized in Supporting Information Table S4. No seroconversion 
was observed in any group at 3 dpi. Antibodies appeared at 7 dpi and all chickens showed 
seroconversion after 21 dpi in the IBV, NDV, and AIV infection groups. The chickens in the 
control group showed negative serum antibody responses against IBV, NDV, and AIV. 

IBV, NDV, and AIV infections were verified using real-time RT-PCR. As shown in Fig. 1, 
negative results were obtained following viral detection in the tracheal tissues from chickens 
in the control group. In the IBV ck/CH/LDL/97I-infected group, the IBV genome was 
detected at 3, 7, 12, and 21 dpi. The viral load gradually decreased during the experimental 
period after the viral copy number peaked at 3 dpi (Fig. 1 A). Similar results were obtained 
from group H120, except no viral genome was detected at 21 dpi (Fig. IB). In the NDV La 
Sota- and AIV H9-infected groups, viruses were detected at 3 dpi and viral copy number 
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peaked at 7 dpi, whereas the results at 12 and 21 dpi were negative (Fig. 1C and D). The 
differences in viral kinetics among IBV ck/CH/LDL/97I, H120, NDV La Sota and AIV H9 
were demonstrated in Fig. IE. Taken together, these results confirmed the successful virus 
infections to SPF chicken. 

3.2 2D-DIGE profiles of tracheal tissues following infections with avian respiratory 
viruses (IBV, NDV, and AIV H9) 

To obtain a detailed comparison of the differential protein expression profiles, the 
tracheal proteins from avian respiratory virus (IBV ck/CH/LDL/97I, HI20, NDV La Sota, 
and AIV H9)-infected and mock-inoculated chickens were extracted for quantitative and 
comparative analysis via 2D-DIGE at 3, 7, 12, and 21 dpi. Comparative analysis revealed that 
total 53,44, 70, and 33 spots were significantly differentially expressed at 3, 7, 12, and 21 dpi 
respectively, which are indicated by arrows in the corresponding representative 2D-DIGE 
map (Fig.2A, B, C and D). In the IBV ck/CH/LDL/971-, H120-, NDV La Sota-, and AIV H9- 
infected groups, and between ck/CH/LDL/971- and H120-infected groups, 72, 65, 63, 86 and 
63 differentially expressed protein spots were identified, respectively (Table 1 and Supporting 
Information Figure SI). In general, these data showed that the expression of more host 
protein spots were changed in the chicken tracheal tissues after virulent IBV ck/CH/LDL/971 
and AIV H9 infections, comparing with those of the vaccine strains IBV HI 20 and NDV La 
Sota strains. 
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3.3 Identification and classification of differentially expressed proteins 

To identify the differentially expressed proteins in viral infected tracheal tissues, 200 
spots with thresholds >1.2-fold (Supporting Information Table S2) were excised from 
preparative gels and subjected to MALDI-TOF/TOF MS identification. Detailed information 
of the identified proteins is provided in Supporting Information Table S5 and Supporting 
Information Fig. S2. Totally, 200 differentially expressed protein spots corresponding to 75 
altered cellular proteins were successfully identified; the biological functions of the 
differentially expressed proteins were summarized in Supporting Information Table S6. 
Number of differentially expressed spots and proteins induced by IBV ck/CH/LDL/97I, H120, 
NDV La Sota, and AIV H9 infections were summarized in Table 1. Among these identified 
proteins, some were identified from two or more separated spots, indicating that they may be 
differently modified ( e.g . phosphorylation of ANXA1 and HSPB1) in the process of host 
response to IBV, NDV and AIV H9 infections (Fig. 2). In order to verify the expression of 
phosphorylated HSPB1, tracheal proteins from sacrificed birds of each group were examined 
by Western blot analysis with specific antibodies to phospho-HSPBl (pSer82). The up- 
regulation of phosphorylated HSPB1 in IBV-, NDV- and AIV H9-infected groups at 7dpi and 
12dpi were indentified (Fig. 3A and 3B); these results were consistent with the expression 
pattern of HSPB1 spots showed by 2D-DIGE, confirming the expression of phosphorylated 
HSPB 1 during IBV, NDV and AIV H9 infections. 

To better understand the potential role of host factors involved in viral infection, 
differentially expressed proteins were classified based on gene ontology annotations derived 
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from the AgBase database ( www.agbase.msstate.edu/ ). As illustrated in Fig. 4A, the biological 
process of proteins differentially expressed in IBV ck/CH/LDL/97I-, H120-, NDV La Sota-, 
and AIV H9- infected groups and between IBV ck/CH/LDL/97I- and H120-infected groups 
(Fig. 4B) all mainly involved in anatomical structure development, biosynthetic process, 
transport, signal transduction, protein complex assembly, cell death, cell differentiation, 
translation, cellular protein modification process, response to stress, cytoskeleton 
organization, growth, macromolecular complex assembly, immune system process, vesicle- 
mediated transport, embryo development, cell proliferation, transmembrane transport and 
various metabolic associated process. These analyses indicated that all of the four viruses 
affect host global response similarly. 

In order to highlight the differences of proteome changes in trachea induced by different 
viruses, a side-by-side comparison of the differentially expressed proteins in trachea between 
different groups was performed (Table 2 and Supporting Information Table S7). The 
comparisons revealed that alteration of certain proteins and corresponding biological process 
appeared to be general or specific to host response in the process of different virus infections. 
For example, EN01, EN03, PGAM1, ANXA1, CKB, GST, HSPB1, HBAA, HBAD, LMNA, 
akr and PRDX1 were commonly altered in all viruses incubated groups, likely suggesting that 
these proteins were involved in and had similar functions in response to the infection of the 
three avian respiratory viruses. GAPDH, LDHA, PPIA, PRDX6 and HADH were altered in 
pathogenic IBV ck/CH/LDL/97I- and AIV H9- infected groups, while PSMB2 was more 
abundant in vaccine strains IBV HI 20- and NDV La Sota- incubated groups. And we also 
found that some proteins were specifically affected by IBV ck/CH/LDL/97I (S100A11, 
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ANXA8 and PHB), H120 (UQCRFS1), NDV La Sota (rbf and VDAC1) and AIV H9 (RPSA, 
CTHRC1, PSMA6 and TAGLN) infections, respectively (Table 2). 

Expectedly, some proteins identified in this study have also been reported in proteomic 
analyses of cells or tissues infected with IBV or AIV (Table 3). For example, up-regulation of 
LGALS1, PPIA, PRDX4 and HSPB 1 has been found in AIV H9- infected human cells [11]. 
ANXA2, LGALS1, EIF5A2 and PPIA were also up-regulated in IBV-infected avian DF-1 
cells [12]. And down-regulation of ANXA1 was reported in IBV-infected chicken trachea in 
our previous works [14, 15]. The common alteration of these proteins in different models of 
infection with the same virus of different strains suggests that these proteins certainly 
involves in the course of corresponding virus infection. 

3.4 Transcriptional profiles of differentially expressed proteins 

Transcriptional alteration analysis of ten selected genes at 3, 7, 12, and 21 dpi were 
achieved by quantifying the mRNA transcripts, in which the housekeeping gene 28S rRNA 
was used as an internal control (Fig. 5A). In this study, primers designed for transcriptional 
alteration analysis were located at positions across intron-exon boundaries to avoid the 
interference of genomic DNA, so quantitative comparison of total amount of mRNA 
transcripts of each gene was conducted. In general, the trends in the change in mRNA 
abundance of selected genes at four time points (Fig. 5A) were similar to the change patterns 
of their corresponding proteins showed by 2D-DIGE (Fig. 5B and Supporting Information 
Table S2). And degrees of change in mRNA abundance of ANXA1, HSPB1, S0D1, EIF5A2, 
LGALS1 and PRDX1 were more obvious than the abundance change of corresponding 
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proteins showed by 2D-DIGE. This is not unexpected because post-transcriptional 
mechanisms, including protein translation, post-translational modification, and degradation 
may also contribute to the abundance change of these proteins during virus infection [28]. 
These data provided transcriptional information complementary to the differentially 
expressed proteins revealed by proteomics analysis. 

3.5 Validation of the differentially expressed proteins by western blot analysis 

To confirm the dynamic alterations of protein abundance during ck/CH/LDL/97I, HI 20, 
La Sota, and H9 infections, four protein candidates ANXA1, PRDX1, ANXA2 and HSPB1 
were chosen for western blot analysis with specific antibodies (Fig. 3A and 3B). The 
abundance changes of ANXA1 and PRDX1 (both altered at all four time points); ANXA2 
(altered only at 12dpi) and HSPB1 (altered at 3, 7, and 12dpi) were consistent with the 
expression changes showed by 2D-DIGE (Fig. 5B and Supporting Information Table S2). 
Altogether, Western blot validated the results obtained from proteomics analyses during IBV, 
NDV and AIV H9 infections. 
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4 Discussion 

In this study, host responses against different avian respiratory viruses which result in 
distinct disease phenotypes were analyzed using comparative proteomics. IBV 
ck/CH/LDL/97I and AIV H9 are two kinds of pathogenic viruses for chicken causing 
significantly clinical signs, while, IBV HI20 and NDV La Sota are avirulent vaccine strains 
inducing effective immune response. Our proteomic analyses of the alterations in protein 
expression induced by different viruses under the same physiological conditions providing 
more specific information of molecular pathogenesis and virus-host interactions of these 
different viruses. Although GO functional annotation and classification showed that all four 
viruses affect host global cellular processes similarly, the trends and degrees of alteration of 
protein expression were not identical. What apparent from the results is that pathogenic IBV 
ck/CH/LDL/97I and AIV H9 induced slightly more changes in protein expression than the 
vaccine strains IBV HI20 and NDV La Sota in chicken trachea. In addition, the differences in 
protein expression changes induced by IBV ck/CH/LDL/97I, HI20, NDV La Sota, and AIV 
H9 indicated that different interactions likely occurred in the host response to infection by 
these viruses and resulted in the different pathogenicity. This is the first report of comparative 
proteomics analysis performed to study the host response to infections of different avian 
respiratory viruses; the information obtained herein might provide critical clues for future 
studies. 
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4.1 Inflammatory responses in IBV, NDV and AIV H9 infections 

IBV, NDV and AIV H9 subtype viruses can initially infect the chicken tracheal 
epithelium and elicit inflammatory responses associated with the disease [16, 29, 30]. In this 
study, several proteins involved in inflammatory responses were identified including ANXA1, 
PPIA, LGALS1 and PRDX4. 

Our previous works have been reported down-regulation of the annexin A1 (ANXA1) in 
IBV-infected trachea [14, 15]. In this study, the endogenous anti-inflammatory protein 
ANXA1 [31] was found down-regulated in both pathogenic IBV ck/CH/LDL/97I- and AIV 
H9-infected groups at 7dpi and in all virus-infected groups at 12dpi, suggesting that the 
decreased expression of ANXA1 may be associated with the inflammatory phenotype during 
the course of avian respiratory viruses infections. More notably, ANXA1 was identified in 
more than one spots in our study, suggesting that ANXA1 phosphorylation during IBV, NDV 
and AIV H9 infections may associates with their pathogenesis [32]. The peptidyl-prolyl 
isomerase (PPIA) has been shown to elicit inflammatory responses and inhibit CD4 + T cell 
activation through its interaction with Itk [33]. We noticed that pro-inflammatory PPIA was 
only up-regulated in pathogenic IBV CK/CH/LDL/97I- and AIV H9-infected groups; there 
was no alteration of PPIA expression in trachea of chickens infected with IBV HI 20 and 
NDV La Sota. This difference may indicate that pathogenic IBV ck/CH/LDL/97I and AIV H9 
trigger more severe inflammatory responses and immune suppression than IBV HI 20 and 
NDV La Sota vaccine strains, which associated with their different pathogenetic processes. In 
addition, PPIA can facilitate replication of human immunodeficiency virus-1 (HIV-1) [34] 
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and coronaviruses including IBV [35]. However, PPIAcan inhibit the influenza A viral 
replication through enhancing the degradation of Ml protein [36]. Chicken PPIA also showed 
anti-influenza virus activity [37]. Up-regulated PPIA has been detected in human cells 
infected with AIV H9N2 [11]. Hence, we speculate that the up-regulation of PPIA may 
contribute to the replication and pathogenesis of IBV ck/CH/LDL/97I and participate in host 
defense against AIV H9. 

Peroxiredoxin 4 (PRDX 4) can protect pancreatic beta-cells against injury by 
suppressing oxidative stress and inflammatory signaling activation [38]. Galectin-1 (LGALS1) 
is an important contributor to immune homeostasis and exhibits anti-inflammatory activity 
during infection [39]. We found that anti-inflammatory PRDX4 and LGALS1 were generally 
up-regulated in IBV ck/CH/LDL/97I-, NDV La Sota- and AIV H9-infected groups at 12dpi. 

In the case of avian respiratory viral infections, the inflammatory process diminished at the 
recovery stage [16]. Hence, the up-regulation of anti-inflammatory PRDX 4 and LGALS1 
may be involved in tracheal epithelium recovery from inflammation damage induced by IBV 
ck/CH/LDL/97I, NDV La Sota and AIV H9. Besides as an immunomodulator, galectin-1 can 
function as pattern recognition receptor (PRR) to facilitate pathogen recognition [39]. The 
interaction of galectin-1 with viral envelope glycoproteins can inhibit infectivity of AIV and 
Nipah virus [40, 41]. Galectin-1 promotes HIV-1 infection through stabilization of viral 
adsorption [42]. Moreover, up-regulation of galectin-1 was also found in human cells infected 
with AIV H9N2 [11], suggesting that galectin-1 may participate in innate immune response to 
AIV H9 infection. Our results also indicate that galectin-1 potentially involved in IBV 
ck/CH/LDL/97I and NDV La Sota infections, whether galectin-1 binds to the surface of IBV 
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and NDV, and functions as effective PRR or promotes infection by facilitating viral 
adsorption need to be further elucidated. 

4.2 IBV, NDV and AIV H9 may induce alterations of nuclear structure and hijack 
host translational machinery 

Nucleus associated proteins, including nuclear mitotic apparatus protein (NUMA), 
lamin-A (LMNA), and eukaryotic initiation factor 5A-l,-2 (eIF5Al, 2), were identified in the 
present study. Though the primary site of replication of IBV, NDV and AIV is the cytoplasm, 
IBV nucleocapsid protein, NDV matrix protein, and AIV nucleoprotein can localize to the 
nucleolus [43-45]. In this study, LMNA was universally down-regulated in all virus-infected 
groups and NUMA was down-regulated in IBV ck/CH/LDL/97I, HI20- and AIV H9-infected 
groups. NUMA and LMNA are components of nuclear matrix and nuclear lamina, down- 
regulations of NUMA and LMNA might implicate that IBV, NDV and AIV H9 infections 
commonly result in disruption of nuclear architecture due to nuclear entry and accumulation 
of viral proteins in the nucleolus. Down-regulation of LMNA has also been detected in 
influenza virus-infected cells [10] and IBV-infected trachea [15]. 

As well known that virus may hijack host translational machinery by regulating distinct 
elFs [46-49]. Previous reports showed eIF5 A can regulate translation and replication of HIV- 
1 [48], and eIF5A mutants inhibit HIV-1 replication [49]. In this study, eIF5Al was down- 
regulated in IBV HI 20-, NDV La Sota- and AIV H9-infected groups, while, eIF5A2 was up- 
regulated in IBV ck/CH/LDL/97I-, NDV La Sota- and AIV H9-infected groups. We speculate 
that IBV, NDV and AIV H9 may interfere with the host protein synthesis by down-regulating 
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constitutively expressed eIF5Al. While the up-regulation of eIF5A2 facilitates viral mRNAs 
translation in infected cells; the expression of eIF5A2 is extremely low in normal conditions. 
Furthermore, the differences in expression of eIF5Al and eIF5A2 indicate that the two IBV 
strains differentially hijack host translational apparatus for their own benefits. Our results 
provide additional evidences that RNA viruses can interact with the nucleolus to hijack 
cellular functions and recruit nucleolus proteins to facilitate virus replication. 

4.3 IBV, NDV and AIV H9 manipulate host UPP through different strategies 

The ubiquitin-proteasome pathway (UPP), a cellular ubiquitin-dependent protein 
degradation system, has been implicated in viral infection cycle and virus-host interactions 
[50]. Inhibition of the UPP has been reported to affect the replication of several viruses, 
including influenza virus [51] and coronavirus [52]. Moreover, the interaction of HIV-1 Tat 
protein with proteasomal [3 subunits contributes to the immune escape of HIV-1 [53]. 
However, the potential roles of host UPP during IBV and NDV infections and precise 
interaction mechanisms between UPP and AIV are not well elucidated. Previous reports 
showed that proteasome subunit alpha and beta types were up-regulated in IBV-infected DF-1 
cells [12]. In this study, proteasomal subunit alpha type-6 (PSMA6) was up-regulated only in 
AIV H9-infected group, whereas proteasomal subunit beta type-2 (PSMB2) was commonly 
up-regulated in IBV HI 20- and NDV La Sota-infected groups and more abundant in H120- 
infected group than ck/CH/LDL/97I-infected group, indicating that different viruses (IBV, 
NDV, and AIV) and different strains of the same virus (ck/CH/LDL/97I and HI20) may 
evolve different strategies to utilize UPP to facilitate replication and infection by specific 
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interaction with distinct proteasomal subunits. 

4.4 Different roles of ANXA2 in IBV, NDV and AIV H9 infections 

Annexin A2 (ANXA2) has been shown involved in exocytosis, cell proliferation, 
membrane trafficking and apoptosis [54], the organized ANXA2 membrane domains are 
necessary for virion budding from infected cells [55]. Notably, many viruses utilize ANXA2 
for effective infection and replication, such as respiratory syncytial virus (RSV) [56], HIV 
[57] and HCV [58]. In this study, ANXA2 was commonly up-regulated during IBV HI 20, 
NDV La Sota and AIV H9 infections. We speculate that the up-regulated ANXA2 may 
support replication of AIV H9 by incorporated into influenza virus particles to enhance 
influenza virus replication via mediated activation of plasminogen [59]. However, up- 
regulated ANXA2 might act as an anti-viral protein to inhibit replication of IBV HI 20 as it 
was reported that ANXA2 reduces viral frameshifting efficiency and successful replication of 
IBV as a cellular RNA-binding protein [54], Up-regulation of ANXA2 has also been found in 
avirulent IBV-infected DF-1 cells [12]. In our previous work, the abundance of ANXA2 was 
only increased following attenuated IBV ck/CH/LDL/97I Pus infection, no significantly 
abundance changes of ANXA2 were observed during virulent IBV ck/CH/LDL/97I P 5 
infection, suggesting that the ANXA2 may be specific host response to avirulent IBV strains 
(ck/CH/LDL/97I Pi 15 and H120). The function of ANXA2 in life cycle of NDV needs further 
investigation. 
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4.5 IBV, NDV and AIV H9 infections and stress response proteins 

Heat shock protein 27 (HSPBl) is a multifunctional protein involves in various cellular 
processes and viral infection [60]. Phosphorylated HSPBl has been shown to participate in 
the replication of adenovirus [60], HCV [61], and herpes simplex virus type 1 [62]. In this 
study, the up-regulation of phosphorylated HSPBl during IBV, NDV and AIV H9 infections 
was confirmed by western bolt (Fig. 3). It was reported that RSV induced phosphorylation of 
HSPBl in human bronchial epithelial cells ultimately results in epithelial membrane 
permeability [63]. Avian respiratory viruses infections can result in tracheal epithelium lesion, 
but the severity and persistence of the lesions differed among different viruses [16-18]. 
Therefore, phosphorylated HSPBl which mainly up-regulated in IBV HI20- and NDV La 
Sota-infected groups at 7dpi and in IBV ck/CH/LDL/97I- and AIV H9- infected groups at 
12 dpi may be responsible for the tracheal epithelium lesions at different stages of avian 
respiratory viruses infections. Further studies should focus on the exact roles of 
phosphorylated HSPBl in the replication and pathogenesis of IBV, NDV and AIV H9. 

Viral infections are capable of induce oxidative stress associated with the activation of 
phagocytosis and release of reactive oxygen species which positively modulate immune 
activation, eliminate viral infection and immune-induced cellular injury [64], In this study, 
antioxidant stress proteins including peroxiredoxin (PRDX)-l, 3,4, 6, glutathione S- 
transferase (GST) and superoxide dismutase 1 (SOD1) (Table 2 and Supporting Information 
Table S2) were universally altered in all viruses-infected groups. Most of altered proteins 
were up-regulated, which may be related to relief of cellular oxidative stress induced by IBV, 
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NDV and AIV H9 infections. Moreover, more extensive changes were observed in AIV H9- 
infected group, suggesting that AIV H9 may activate more intense oxidative stress. Among 
these antioxidant stress proteins, PRDX1 and PRDX4 have been reported to protect RSV- 
induced oxidative damage to the nuclear cytoskeletal proteins in epithelial cells [65], while 
down-regulation of antioxidant enzymes (SOD, GST and PRDX) induced by RSV 
contributes to RSV pathogenesis [66]. In addition, PRDX1 can activate secretion of 
proinflammatory cytokines by interacting with TLR4 [67]. PRDX4 plays critical protective 
roles in inflammatory diseases [38]. Similarly, we considered that the up-regulation of 
PRDX1 and PRDX4 during IBV, NDV, and AIV H9 infections may activate host antiviral 
immunity to eliminate viral infection and provide protective cellular mechanisms to reduce 
oxidative stress and inflammatory damage. However, the potential roles of down-regulated 
oxidative stress-associated proteins in the pathogenesis of IBV, NDV, and AIV needed to be 
investigated further. 

In this study, trachea was chosen for global analysis of the host response because it is the 
primary target organ of avian respiratory viruses. The comparative analysis of changes in 
tracheal protein profile induced by IBV, NDV, and AIV H9 in vivo provides insight into 
mechanisms of global host response to different avian respiratory virus infections. However, 
the sites of IBV, NDV, and AIV infection are actually not completely consistent, IBV 
replicates in the epithelium of lower respiratory tract of chicken [68], NDV and AIV H9 
subtype induce lesions in the upper respiratory tract [18, 69]. The roles of different sites of 
respiratory tract in this study had not been effectively measured. Moreover, infections of IBV, 
NDV and AIV H9 subtype viruses can elicit inflammatory responses characterized by 
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recruitment of myeloid or lymphoid cells. Consequently, the cellular types of tracheal tissues 
harvested from mock and virus-infected chickens for proteomics analysis may differ due to 
immune cell infiltration. As a result, the changes revealed by proteomic analysis of some 
proteins such as LGALS1, PPIA and PRDX1 that could also be expressed in immune cells in 
the course of viral infection may be caused or partly caused by cellular heterogeneity 
occurred during the inflammatory response to viral infection (including infiltration by innate 
and adaptive immune cells). This needed to be further investigated. In addition, though a 
number of differentially expressed proteins that have potential roles in revealing the common 
or specific aspects of viral pathogenesis and virus-host interactions were identified, the 
precise roles of these proteins during viral infection had not proved in this study. 
Nevertheless, the similar and different protein expression profiles induced by different avian 
respiratory viruses after infection might be indicators of the essential response of infected 
target cells, provide insights into viral pathogenesis and host antiviral defense, and supply 
useful clues for the development of novel prevention or therapeutics strategies against these 
viruses. 
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Figure 1. The viral load was quantified using real-time RT-PCR. The viral loads in tracheal 
samples of three sacrificial birds from each group at each time point were determined. The 
average viral copy number per 1 pg RNA of each group was calculated through standard 
curve based on ten-fold dilution series of standard template with known concentration. The 
limit of detection of the real-time RT-PCR reaction was 10 2 , 10 1 and 10 1 copies of template 
for IBV, NDV and AIV, respectively. (A) Kinetics of viral loads of IBV ck/CH/LDL/97I. (B) 
Kinetics of viral loads of IBV HI 20. (C) Kinetics of viral loads of NDV La Sota. (D) Kinetics 
of viral loads of and AIV H9. (E) The differences in viral kinetics among IBV 
ck/CH/LDL/97I-, H120-, NDV La Sota-, and AIV H9-infected groups. Data represent means 
of three biological replicates per group. Error bars indicate standard error of the mean, and 
dpi represent days post-inoculation. The symbol * indicate significant differences (* represent 
p <0.05 and ** represent p <0.01) between different groups. 
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Figure 1. 
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Figure 2. The representative gel images showing 2D-DIGE maps labeled the relative 
locations of spots that displayed significant quantitative differential expression (p<0.05 and 
average ratio>1.2 or <-1.2) during different sampling times. (A) 3 dpi; (B) 7 dpi; (C) 12 dpi; 
and (D) 21 dpi. Representative 2D-DIGE images of IB V ck/CH/LDL/97I-, H120-, NDVLa 
Sota-, and AIV H9-infected groups with their respective sampling points were represented as 
Supporting Information Figure SI. Fluor-labeled protein samples were separated on 24-cm 
pH 3-10 linear IPG strips, subjected to SDS-PAGE, and then gels were scanned on a 
Typhoon 9400 scanner, and image analysis was performed with Ettan™ DeCyder Software 
version v6.5. The identified protein spots were marked with arrows and labeled with the 
respective Match ID listed in Supporting Information Table S2. The green and red arrows and 
ellipses indicating HSPB1 and ANXA1 identified from separated spots. 

Figure 2. 



This article is protected by copyright. All rights reserved 



Accepted Article 


www.proteomics-journal.com 


Page 42 


Proteomics 


Figure 3. Western blot confirmation of representative proteins in IBV ck/CH/LDL/97I-, 
H120-, NDV La Sota-, and AIV H9-infected chickens. (A) The immunoblot analysis of 
ANXA1, ANXA2, PRDX1, HSPB1 and phosphorylated HSPB1. (B) The relative 
densitometric intensity of ANXA1, ANXA2, PRDX1, HSPB1 and phosphorylated HSPB1 in 
immunoblot analysis. [3-Tubulin was used as an internal control to normalize the quantitative 
data. 
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Figure 4. Gene Ontology annotation analyses of differentially expressed proteins according 
to their biological process categories. This classification was produced by using the 
GOSIimViewer tool of the Agbase database (http://www.agbase.msstate.edu/) . (A) Biological 
process classification of all indentified proteins. (B) Comparison of biological process 
classification of differentially expressed proteins in IBV ck/CH/LDL/97I-, H120-, NDV La 
Sota- and AIV H9-infected groups and between IBV ck/CH/LDL/97I- and H120-infected 
groups. 

Figure 4. 
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Figure 5. Transcript alteration of selected differentially expressed proteins in IBV 
ck/CH/LDL/97I-, H120-, NDV La Sota-, and AIV H9- infected chicken. (A) Comparison of 
mRNA levels of selected proteins differentially expressed at 3, 7, 12, and 21 dpi. Total RNA 
extracts were prepared from the tracheal tissues of sacrificial birds in control and all virus- 
infected groups at four time points (3, 7, 12, and 21 dpi), measured by real-time RT-PCR. 
Data represent means of three biological replicates per group. Error bars indicate standard 
error. Gene names refer to Supporting Information Table S6. Samples were normalized with 
the expression of the 28S ribosomal RNA gene in the control group. (B) Average ratios of 
selected proteins determined by 2D-DIGE using DeCyder Software. 
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Table 1. Summary of differentially expressed spots and proteins in IBV ck/CH/LDL/97I-, 
H120-, NDV La Sota- and AIV H9-infected groups and between IBV ck/CH/LDL/97I- and 
H120-infected groups 


IBV IBV NDV La AIV IBV 

ck/CH/LDL/9717 HI 20/Control Sota/Control H9/Control ck/CH/LDL/971/ 

Control HI 20 
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proteins 
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a) Total number of differentially expressed spots and protein. 

b) Number of up-regulated spots and proteins. 

c) Number of down-regulated spots and proteins. 
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carbohydrate metabolic process 7 
generation of precursor metabolites and 
energy 2 
translation 1 
protein folding 1 

cellular protein modification process 2 

lipid metabolic process 2 

response to stress 1 

cytoskeleton organization 1 

cell death 3 

catabolic process 3 

biosynthetic process 2 

small molecule metabolic process 3 

membrane organization 1 

transport 2 

vesicle-mediated transport 1 
cell differentiation 2 
growth 1 

anatomical structure development 4 
protein maturation 1 
membrane organization 1 
catabolic process 2 


C0L2A1 cell adhesion 1 

cell death 1 

embryo development 3 

cell differentiation 1 

extracellular matrix organization 1 

homeostatic process 1 

anatomical structure formation involved in 

morphogenesis 1 

anatomical structure development 18 
neurological system process 2 




www.proteomics-journal.com 


Page 49 


Proteomics 




IBV ck/CH/LDLV97PControl 
IBV HI 20/Control 
IBV ck/CH/LDLV97PControl 
NDV La Sota/Control 



IBV HI 20/Control 
AIV H9/Control 




IBV HI 20/Control 
IBV ck/CH/LDL/97I/ H120 
NDV La Sota/Control 
IBV ck/CH/LDL/97I/ H120 



1 

3 


3 


1 

1 


Q* 

0 > 

CJ_ 

This article is protected by copyright. All rights reserved 


THBS1 


aldolase C 
C4 
OGN 
HAPLN1 
NME1 
MY0Z1 


ECU 

GNB2L1 


cell adhesion 1 

cell death 1 

cell proliferation 1 

biosynthetic process 5 

cell differentiation 1 

cellular nitrogen compound metabolic 

process 6 

small molecule metabolic process 6 
anatomical structure development 2 


protein complex assembly 1 

cellular protein modification process 2 

transport 1 

cell cycle 1 

signal transduction 4 

cell death 3 

catabolic process 2 

embryo development 1 

vesicle-mediated transport 1 

cellular nitrogen compound metabolic 

process 1 

nucleobase-containing compound catabolic 
process 1 
growth 1 
locomotion 1 

small molecule metabolic process 1 
anatomical structure development 1 
cell motility 1 
cell division 1 




www.proteomics-journal.com 




ffiV ck/CH/LDLV97I/Control 




IBV HI 20/Control 
NDV La Sota/Control 


1 

2 




This article is protected by copyright. All rights reserved 


50 


Proteomics 


COPE 

S0D1 

SRL 


S100A11 

ANXA8 

PHB 


UQCRFS1 

rbf 

VDAC1 


RPSA 

CTHRC1 

PSMA6 


transport 3 
catabolic process 1 
vesicle-mediated transport 2 
cellular nitrogen compound metabolic 
process 1 

nucleobase-containing compound catabolic 
process 1 

small molecule metabolic process 1 

DNA metabolic process 1 

cellular protein modification process 1 

signal transduction 3 

cell proliferation 1 

biosynthetic process 5 

cellular nitrogen compound metabolic 

process 4 

growth 1 

chromosome organization 1 

generation of precursor metabolites and 

energy 1 

transport 6 

response to stress 1 

cell-cell signaling 2 

catabolic process 1 

biosynthetic process 1 

cellular nitrogen compound metabolic 

process 2 

nucleobase-containing compound catabolic 
process 1 

small molecule metabolic process 2 
neurological system process 3 
transmembrane transport 3 


TAGLN 




Accepted Article 

H S w 


www.proteomics-journal.com 


Page 51 


Proteomics 


Groups incubated different viruses (IBV ck/CH/LDL/97I, HI20, NDV La Sota and AIV H9) containing commonly expressed proteins. 

Biological process of commonly expressed proteins in different viruses infected groups. This annotation was obtained using GOSlimViewer tool at Agbase 
database ( http://www.agbase.msstate.edu/ ). - indicated no GO annotation was acquired. 
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Table 3. Proteins identified both in our study and previous reports in proteomes of the infected tissues or cells induced by AIV and IBV infections* 
*NDV was not included in this table because the proteome changes of cells or tissues induced by NDV have not yet been reported. 
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a) Proteins identified in this study that induced by AIV H9 infection. 

b) Proteins identified both in our study and previous reports induced by AIV infection. 

c) Proteins identified in this study that induced by IBV ck/CH/LDL/97I and H120 infections. 

d) Proteins identified both in our study and previous reports induced by IBV ck/CH/LDL/97I and H120 infections. 
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